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Primordial magnetic fields lead to non-Gaussian signals in the cosmic microwave background (CMB) even at 
the lowest order, as magnetic stresses and the temperature anisotropy they induce depend quadratically on the 
magnetic field. In contrast, CMB non-Gaussianity due to inflationary scalar perturbations arises only as a higher 
order effect. Apart from a compensated scalar mode, stochastic primordial magnetic fields also produce scalar 
anisotropic stress that remains uncompensated till neutrino decoupling. This gives rise to an adiabatic-like scalar 
perturbation mode that evolves passively thereafter (called the passive mode). We compute the CMB reduced 
bispectrum (ftj^i ) induced by this passive mode, sourced via the Sachs-Wolfe effect, on large angular scales. 
For any configuration of bispectrum, taking a partial sum over mode-coupling terms, we find a typical value of 
/i(/i+l)Z3(i3 + l)&j j l ~ 6 — 9 x 10~ 16 , for a magnetic field of Bo ~ 3 nG, assuming a nearly scale-invariant 
magnetic spectrum . We also evaluate, in full, the bispectrum for the squeezed collinear configuration over all 
angular mode-coupling terms and find h(h + l)h(ls + l)&! 1 i 9 ; 3 ~ —1.4 x 10~ 16 . These values are more 
than ~ 10 6 times larger than the previously calculated magnetic compensated scalar mode CMB bispectrum. 
Observational limits on the bispectrum from WMAP7 data allow us to set upper limits of Bo ~ 2 nG on the 
present value of the cosmic magnetic field of primordial origin. This is over 10 times more stringent than earlier 
limits on Bo based on the compensated mode bispectrum. 

PACS numbers: 98.62.En, 98.70.Vc, 98.80.Cq 



I. INTRODUCTION 

The origin and evolution of large-scale magnetic fields in 
the universe is not yet clearly understood. There is observa- 
tional evidence for micro-Gauss strength magnetic fields, or- 
dered on kiloparsec to 10 kiloparsec scales, in galaxies and 
clusters of galaxies yJ-Ql]. There is also evidence for galactic 
scale fields at higher redshift 10,0]. There are only tentative 
indications, to date, of fields not associated with individual 
galaxies or clusters Q8|] . Further, recent 7-ray observations us- 
ing Fermi/LAT data provide hints of a detection or a lower 
bound [10] of B ~ 10~ 16 G on intergalactic scales. Con- 
straints on purely cosmological magnetic fields have also been 
derived from the CMB fiu430ll . big bang nucleosynthesis and 
polarized radiation from extragalactic radio sources (for a re- 
view see e.g. (HH]). 

It is generally accepted that the observed fields require a 
seed magnetic field, possibly of primordial origin, which are 
then amplified by astrophysical processes JH [j, l3l"l - [34ll . A 
strong enough seed may mainly require amplification due to 
flux-freezing which arises during the collapse to form struc- 
tures D3l. On the other hand a weak seed would require con- 
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siderable dynamo action as well ll3il4l l32l - [34ll . As yet no com- 
pelling theory exists for the origin of strong primordial fields. 
Equally, dynamo theories are not without their own difficul- 
ties, in particular, how large-scale dynamos lead to coherent 
enough fields on saturation Thus it is important to keep 
an open mind on the issue of the origin of magnetic fields. 

Large-scale primordial fields could have resulted from 
phase transitions in the early universe, for example during 
the inflationary era I35l - l39ll . An important way to constrain 
or detect such fields is via their imprints on the Cosmic Mi- 
crowave Background (CMB) anisotropies. Considerable work 
has already been done to calculate the magnetic field sig- 
nals in the power spectrum (both temperature and polariza- 
tion) of the CMB anisotropies IHrMEIQ. The possibility 
of non-Gaussian signals in the CMB temperature anisotropies 
has drawn increased attention of late. This possibility is par- 
ticularly relevant when one considers magnetically induced 
anisotropies in the CMB, for the following reason: In stan- 
dard inflationary models, the Gaussian statistics of the quan- 
tum fluctuations in the inflaton field leads to Gaussian statis- 
tics for the CMB temperature field. In such models, any non- 
Gaussian signal in the t emp erature field arises generically due 
to higher order effects ll46l - l5lll . 

The situation is different in the case of magnetically in- 
duced signals. Magnetic stresses depend quadratically on 
the field. Hence, even if the fields are assumed to be Gaus- 
sian, the corresponding stresses are not. This implies that the 
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CMB temperature anisotropy, which these stresses induce, do 
not have Gaussian statistics. Unlike the case of CMB non- 
Gaussianity from inflationary perturbations, in the case of 
magnetic fields, the CMB non-Gaussianity is a leading order 
effect. Thus the study of CMB non-Gaussianity has a special 
significance in the context of probing and detecting primordial 
cosmological magnetic fields. 

Preliminary studies of such non-Gaussian signals in the 
CMB induced by primordial magnetic fields have begun ll52l - 
54], based on earlier calculations of non-Gaussianity in the 
magnetic stress energy B55I1 . These have been limited to the 
scalar mode signals on large scales and further restricted to 
a component referred to as the compensated scalar mode by 
ll40ll . However, the work of Shaw and Lewis i40ll has revealed 
that much larger CMB anisotropies can result from scalar per- 
turbations sourced by the magnetic anisotropic stresses prior 
to neutrino decoupling (also see II56|-d8I1 ). This mode is re- 
ferred to as the passive scalar mode 1 40] . In this paper we cal- 
culate the contribution to the CMB bispectrum from the pas- 
sive scalar mode. As we show here, this contribution greatly 
dominates over the contribution calculated earlier by two of 
us |p2ll (hereafter referred to as SS09), and also calculated by 
113311 - The current work allows tighter limits to be placed on 
primordial magnetic fields. 

In the next section we describe the properties of the pri- 
mordial stochastic magnetic field and in section 3 we describe 
the perturbation induced due to the passive scalar mode. The 
magnetic CMB temperature anisotropy is described in Section 
4 and in its subsections we calculate the three-point correla- 
tion of the scalar anisotropic stress in two different configura- 
tions. In section 5 we calculate the passive mode reduced bis- 
pectrum for the different configurations and in section 6 we 
use the reduced magnetic bispectra to put constraints on the 
primordial magnetic field. Our conclusions are summarized 
in section 7. 



II. PRIMORDIAL STOCHASTIC MAGNETIC FIELD 

We consider a stochastic magnetic field B which is a Gaus- 
sian random field characterized and completely specified by 
its energy power spectrum M(k). In addition we assume the 
magnetic field is non-helical. 

On galactic and larger scales, any velocity induced by 
Lorentz forces is generally so small that it does not lead to 
appreciable distortion of the initial field J59[|60ll . Hence, the 
magnetic field simply redshifts away as B(x, t) = t>o(x)/a 2 , 
where, bo is the magnetic field at the present epoch (i.e. at 
z — or a — 1). 

We define b(k) as the Fourier transform of the magnetic 
field bo(x). The energy power spectrum is defined by the 
relation (6,(k)6*(q)) = (2n) 3 S(k - q)P y (k)M(fc), where 
fjj(k) — (5ij — kikj/k 2 ) is the projection operator ensur- 
ing V • b = 0. This gives (bg) = 2 J(dk/k)Al(k), where 
A^(fc) = k 3 M (k) / (2tt 2 ) is the power per logarithmic inter- 
val in fc-space residing in the stochastic magnetic field. 

As in [61], we assume a power-law magnetic energy spec- 
tra, M(k) = Ak n that has a cutoff at k = k c , where k c is the 



Alfven-wave damping length-scale lf59l l60ll . We fix A by de- 
manding that the variance of the magnetic field smoothed over 
a 'galactic' scale, kc — lftMpc" 1 , (using a sharp fc-space fil- 
ter) is Bq. This gives, (for n > —3 and for k < fc c ) 

The magnetic spectral index is restricted to —3 < n < 
—3/2. Blue spectral indices are strongly disfavoured by many 
observations like the CMB power spectra 12811 . 112911 and espe- 
cially by the gravitational wave limits of [62]. 

III. PASSIVE SCALAR MODE 

The stress tensor (space-space part of the energy- 
momentum tensor) for magnetic fields in terms of the present 
day magnetic field value bo is 

^ (x) = 4iG^ (x) ^"^ (x) ^ (x) ) (2) 

In Fourier space, the product of magnetic fields becomes a 
convolution 

Sj(k) = J & l (q)Mk-q)d 3 q (3) 

^ (k)= 4iG^ (k) ^~^ (k) )- (4) 

This can be expressed in terms of the magnetic perturbations 
to the energy-momentum tensor as 

T;(k)= Pl (A B s} + n B }) (5) 

where A b and Hb } are the perturbations in the energy density 
and anisotropic stress, respectively, as defined in [40] and p 7 
is the radiation pressure. 

The magnetic stresses are non-linear in the field but we as- 
sume that they are always small compared to the total energy 
density and pressure of the photons, baryons etc. Thus allow- 
ing a purely linear treatment of the perturbations. Hence the 
scalar, vector and tensor perturbations decouple and evolve in- 
dependently and here we focus on the anisotropic scalar per- 
turbations Hs(k) which are given by applying the relevant 
projection operator to T ? (k) Il55ll . 

U B (k) - ~ (kik - ^) U% (6) 

Note that U B (k) of fH is equal to -r s (k) of [0. 

Such an anisotropic stress term also arises due to neutri- 
nos after they decouple. The net anisotropic stress tensor acts 
as the source term in the scalar Einstein's equations that lead 
to the Bardeen equation for the potential. Prior to neutrino 
decoupling, the only source of anisotropic stress is the mag- 
netic field. Once the neutrinos decouple, the anisotropic stress 
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due to neutrinos also contributes but with an opposite sign 
to that of the magnetic field, thus compensating the contribu- 
tion from the magnetic field B18||. After compensation there 
are two modes of perturbation [40]. The first one, called the 
passive mode is an adiabatic-like mode but has non-Gaussian 
statistics. It grows logarithmically in amplitude between the 
epochs of magnetic field generation and neutrino decoupling 
but then evolves passively after neutrino decoupling. This be- 
haviour has also been confirmed in Il57ll in the context of deriv- 
ing the magnetic Sachs-Wolfe effect for a causally generated 
primordial magnetic field. The second, more well-studied per- 
turbation Il23l426ll . is called the compensated mode. 

The final curvature perturbation due to the passive mode is 
given by equation (86) of Shaw & Lewis l40ll . 



C = C(tb) 



1 
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(7) 



in the conformal Newtonian gauge. The evolution of the cur- 
vature perturbation has also been discussed (in synchronous 
gauge) in B56I1 for the case of an extra source of anisotropic 
stress cancelling the neutrino anisotropic stress. The role of 
anisotropic stresses on CMB has also been discussed by 16311 . 
Here R u ~ 0.4 is the fractional contribution of neutrino en- 
ergy density towards the total energy density of the relativis- 
tic component. The fractional contribution of radiation energy 
density is i? 7 = 1 — R u ~ 0.6. The £ (tb) term represents 
all primordial contributions to the curvature perturbation. The 
log term shows the growth in the curvature between the epochs 
of magnetic field generation tb and neutrino decoupling t„. 
This log term is always bigger than 10 for different choices of 
tb so we can neglect the ((5/87?^) — 1) term with less than 5% 
error. After neutrino decoupling, the anisotropic stress ILg is 
compensated and there is a remnant adiabatic mode with am- 
plitude given by [40]. 
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This so-called passive mode evolves passively like an adia- 
batic perturbation but the statistics of n B are non-Gaussian 
unlike the standard inflationary adiabatic mode. Also, unlike 
the compensated mode, the passive mode amplitude depends 
logarithmically on the epoch of magnetic field generation. 
Therefore, limits on the primordial magnetic field strength 
arising from this passive mode will be somewhat sensitive to 
the assumed model and epoch of magnetic field generation. 



IV. MAGNETIC CMB ANISOTROPY AND THREE-POINT 
CORRELATION 

We consider the CMB temperature anisotropies that are 
sourced by the scalar passive mode. On large angular scales 
the important contribution is via the magnetically induced 
Sachs-Wolfe effect 



AT 



1 
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(n) = -3>(x -nD*) = -((x -nD*) (9) 



where the second equality is from, for example, Eq. (10.42) 
of [64]. The unit vector n gives the direction of observation 
and D* is the (angular diameter) distance to the surface of last 
scatter. Then employing Eq. ^ we get 



AT 



(n) 



TZ p U B (x Q -nD*). 



where we define 
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In SS09 we had calculated the bispectrum due to the com- 
pensated scalar mode for which AT/T = (i? 7 /20)A B I65L 
Equation (TTOb for the passive scalar mode leads to a AT/T 
which is larger by a factor A = (4/3) In {t v /tb) (Hb/^-b)- 
Assuming that the fractional perturbations A b and Hb are of 
the same order, we get A of order 50 for tb ~ 10 14 GeV. As 
we will see below for the magnetic spectra we consider, this 
alone leads to an enhanced passive scalar mode contribution 
to the bispectrum of order A 3 w 10 5 and thus stronger limits 
on the primordial magnetic field by a factor A ' 5 . We assume 
instantaneous recombination which is a good approximation 
for large angular scales. There could also be additional inte- 
grated Sachs-Wolfe (ISW) contributions to AT/T(n). 

The CMB temperature fluctuation in a direction (n) at the 
observer can be expanded in terms of the spherical harmonics 
to give 



AT(n) ^ 

— — = 2^ aimYimtn). 

Ira 



(12) 



where, for the passive scalar mode, 



4ir(-i)< 



d 3 k 
(2^)3 



TZ p U B (k) i ; (fcT>*)Fj* m (fe). (13) 



Here Hb (k) is the Fourier transform ofUs (x) . 

The non-Gaussianity in the CMB temperature anisotropy 
can be evaluated by calculating its 3-point correlation function 

7TI 7TI 7YI 

(in harmonic space), called the bispectrum, B l I t 2 3 . We 
assume that the magnetic perturbations are uncorrelated with 
the primary (inflationary) perturbations in the CMB. 
In terms of the a; m 's the CMB bispectrum is given by, 



(^limi ^l^ra^ ^l^ms ) ■ 



(14) 



From Eq. (fT3l we can express B™\ ™ 2 ™ 3 



as 



Ift-0' 



.1=1 



d 3 kl 
'2^2 



with f 123 defined as, 

C 123 = <n B (A0n B (AgrL3(fc 3 )>. 



C123 

(15) 
(16) 



The function ( 123 is the three-point correlation of IT B (fe) and 
involves a 6-point correlation function of the magnetic fields. 
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Using Wick's theorem it can be written as a function of mag- 
netic spectra in an analogous manner to [55] and Eq. (7) of 
SS09, C 123 = S(ki + k 2 + fc 3 ) V123. where 



^123 = J d3sM ^ kl + »\)M{a)M(\a - fc 3 |) 

x (Jn B n B n B ). (17) 



where J-ti b ti b ti b is the angular component of the three-point 
correlation f 123 of the scalar anisotropic stress IT^fc) in 
mode-coupling integral V>i23- 



A. Three-Point Correlation of Scalar Anisotropic Stress 

This angular component J r n B n B n B is given by a 58- 
term expression in [55] derived by applying relevant pro- 
jection operators to extract the scalar part of the full bis- 
pectrum i.e the six-point correlation of the magnetic fields. 
In this particular case, the operator is given by Aijkimn — 
(-1) 3 Qij{k 3 )Qki(ki)Qmn(k 2 ) where Q ab (k) = S ab - 
(3/2)P a b(k) and P a b(k) = S a b — k a kb is the projection op- 
erator. We present the full angular component (in our notation 
we have absorbed a factor of 8 multiplying T in 05511 into our 
definition of Jn B n fl n B ): 

6 

with 



■^nBn B n B - 9 
■^ii B n B n B = 

^Ln B n B = +1 2 + P 2 + 9(0? 3 + 0is + ° 2 i2) + 3(«3 + «? + a% + (3 2 3 + ft + f3 2 2 + 7s 2 + if + ll)) 

Fn B n B n B = ^ 3 ^m(^373 + Pili + P2I2 + ^7) + 7 («373 + "i7i + "272) + P(a 3 (3 3 + + a 2 f3 2 ) 
+38 13 {a 3 a 1 + /3 3 /3i + 7371) + i0 2 z{a z a 2 + /3 3 /3 2 + 7372) + W 12 {aia 2 + + I1I2) 

+9013023^12 

■fn B n B n B = 3 (iPaaPa + /3m«i7i + P1P2I2 + 3 (/20i 3 /3 3 7i + 7023«372 + /30 12 ai/3 2 ) 
+3(a 3 /3 3 (ai/3i + a 2 fi 2 ) + ai7i(a 3 7 3 + a 2 7 2 ) + ^272(^373 + Pill)) 
+9(6>i 3 6»237i72 + Si 3 9i 2 hP 2 + 23 0i 2 a 3 ai) 

■?n B n B n B = ^ 9 f/2a 3 /3 3 ai 7 i + ja 3 /3 3 /3 2 j 2 + /3ai7i/3 2 72 + 3(0i 3 /3 3 7i/3 2 72 + 023«3ai7i72 + Si 2 a 3 p 3 aiP 2 ) 
J^ B n B n B = 27a 3 /3 3 ai7i/5272. (19) 



where 

f3 = (s • fe 3 — s), 7 = (s ■ fci + s), /l = (k 3 - s ■ fei + a). 

(20) 

where the hat on a vector denotes its unit vector. 
Also, 

a a = k a ■ s, p a = k a ■ k 3 - s, 7 a = fc a ■ fei + s (21) 
and 

6 ab = k a -k b . (22) 

where our angle definitions are consistent with lf55ll and 
slightly different to SS09. 



For simplicity, we evaluate this Tu B n B n B expression to 
find the bispectrum in the following two cases below. We ex- 
pect that the magnitude of the bispectrum will be of similar 
order for a general case. 

Case I - Consider any bispectrum configuration but include 
only the s-independent terms (constant for any particular con- 
figuration) in the Jn B n B n B expression. Then the s-integral 
in Eq. ( TTTb is performed without reference to any particular 
bispectrum configuration. 

Case II - The squeezed collinear configuration where we 
calculate fully the mode-coupling integral over all angular 
terms in the Jn B n B njj expression. 
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TABLE I. The sum of s-independent terms m = Jn n s n B m f° ur 
different configurations (fei, k 2 , fc 3 ) for evaluating the Case I bis- 
pectrum. 



Configuration 


(fei, fe 2 , fc 3 ) 


(#12, #2 3 , #1 3 ) 


m 


local 
isosceles 


ki ~ k 3 
k 2 < ki, k 3 


(0,0,-1) 





equilateral 


kl ~ k 2 ~ fc 3 


(1 1 1\ 

V 2 ' 2 ' 2 / 


— O.DZO 


squeezed 
collinear 


fci ~ fc 3 
fc 2 < fci , fc 3 
fei = &2 = — fc 3 


(1,-1,-1) 




midpoint 
collinear 


fa ~ fei ~ ^ 

fei = fc2 = -fc 3 


(1,-1,-1) 


-9 



il For the squeezed collinear configuration case, ■7 r n B n B n B P^ks up 
another term p? ~ 1 



B. Case I: Evaluation with s-independent Terms in Jn B n B n B 

The s-independent terms include #12, #23, #13 that are con- 
stant for a given (fei,fe2,fe 3 ) configuration. For each con- 
figuration considered, we calculate the sum of the five s- 
independent terms in Eq. dT9b for Jn B n B n B 



n B n B n f 



-9 + 9 1 



'23 



27 ( 



(23) 

We give in Table I the values of m = J~ji b ii b iib i^ 1 ' ^2, ^3) 
for specific configurations (fci, ki2, k^). Note that m happens 
to vanish exactly for the local isosceles configuration. Hence 
the Case I mode-coupling integral reduces to 



3 3 m 



123 



(47rp 7 ) 3 



, 7 , (47Tp ) 3 
where po is the present day energy density of radiation and 



(24) 



1= J d 3 sM(\k 1 + s\)M(s)M(\s - fc 3 |) 



2,t. V I dn / dss n+2 (s 2 + k\+ 2skxv) 2 
-1 Jo 



x (s 2 + kl - 2sk 3 fi) 



(25) 



where v — k\ ■ s and j« = fc 3 • s. 

We perform the mode-coupling integral using the technique 
discussed in 01 El Bui. As 771 vanishes for the local isosce- 
les configuration, ^123 is zero for this configuration. For 
the equilateral and squeezed collinear configurations we split 
the s integral into two sub-ranges < s < k\ ~ fc 3 and 
s > fci ~ &3 We then approximate the integrands by assum- 
ing s <C fci ~ fc3 and s 3> fc x ^ fc 3 for the respective sub- 
ranges. For the midpoint collinear configuration we split the 
s integral into two sub-ranges < s < fci and s > 2 fci ~ fc 3 



while neglecting the very small contribution from the middle 
sub-range fci < s < 2 fci ~ fc 3 . Again, we approximate 
the integrands by assuming s fci and s 3> 2 fci ~ fc 3 for 
the respective sub-ranges. To derive numerical estimates of 
the bispectrum and magnetic field strengths we will focus on 
nearly scale-invariant spectra (which can be produced by an 
acausal mechanism like inflation), i.e n — > —3, which yield 



4irA J 



K 



2n+3 , 



4ttA 



3 ft l 



n + 3 
2n ~'~ 3 fc n 



'371+3 



3n + 3 



n - 



(26) 



The latter equation is obtained in the limit n — > — 3 where we 
can neglect the terms with (3n + 3) _1 compared to (n + 3) _1 . 
Hence the mode-coupling integral for Case I - taking only s- 
independent angular part - is 



Vl (27) 



Here we have defined Va, the Alfven velocity in the radiation 
dominated era as 160] 



Va = 



Bo 



(ie^o/3) 1 / 2 

with £_9 = (S /10~ 9 Gauss). 



3.8 x lCr 4 £_ 9 , 



(28) 



C. Case II: Squeezed Collinear Configuration - All Angular 
Terms 

For Case II - We take the squeezed collinear bispectrum 
configuration as fci ~ fc 3 and &2 <C fci , fc3 and fci = fc^ = 
— fc 3 . One wavevector (£2) is negligibly small compared to 
the other two which are also in exactly opposite directions 
and of a similar magnitude. This affords considerable simpli- 
fication and reduction of the angular terms given by Eq. ( fT9l ) 
in doing the mode-coupling integral. Using fc 3 ps —fci we see 
that 



/3 = (S ■ fe 3 — s) PS (S ■ —fci — s) PS — (s ■ fci + s) 



-7- 
(29) 



/1 = (fc 3 — S ■ fci + s) ps (—fci — s ■ fci + s) 
ps -(fcT+s -ki+s) ps -1. (30) 



Q'l = Qf 2 



-a 3 



-/'• 



ft = = -03 = P- 



7i = 72 = -73 = 7- 



'12 — J-j "23 — "13 



(31) 



(32) 



(33) 



(34) 
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These relations substituted into Eq. dT~9b reduce the 58 terms 
to a 10-term angular expression for the squeezed collinear 
configuration 



-8 + /3 2 + 9 (fi 2 + 2 7 2 ) + 6^/3 7 + 3/3 2 7 2 

-9 7 2 (3a/ 2 + 7 2 ) - 18/i^7 3 + 27//V- (35) 

We again perform the mode-coupling as discussed above and 
split the s integral into two sub-ranges < s < ki ~ k$ and 
s > ki ~ k%. Once again we approximate the integrands by 
assuming s <C fci ~ £3 and s 3> fci ~ fc 3 for the respective 
sub-ranges to give 



^123 



I 



(4ttp 7 ) 3 



1 = 



(36) 



where 



1= J d 3 sM(\k 1 + s\)M(s)M(\s -fc 3 |) 

x [-8 + /3 2 + 9 (^ 2 + 2 7 2 ) + 6/487 + 3^V 
- 9 7 2 (3a/ 2 + 7 2 ) - 18//^7 3 + 27m 2 7 4 ] 



(37) 



V. PASSIVE SCALAR CMB BISPECTRUM 

Using Eq. (l40l > in Eq. (Tl3T > we can evaluate the CMB bis- 
pectrum for the passive scalar mode for both Case I and II. 
The algebraic steps are the same as those for the compensated 
scalar mode in SS09. We express the delta function present in 
£ 123 in its integral form S (fc) = (l/(27r) 3 ) J d 3 x exp(ifc ■ x), 
use the spherical wave expansion of the exponential terms, 
substitute it into Eq. dl51 l. and integrate over the angular parts 
of (fci, fc2, fc3, x). This algebra is also very similar to that in 
the calculation of the primordial bispectrum [49]. Then it be- 

T71> TTL TYL 

comes possible to write the bispectrum B x \ 1 2 3 , in terms of 
a reduced bispectrum bi 1 1 (also referred to as the Komatsu- 
Spergel estimator 114711 ) as 



(42) 



where 



3 

J] / kfdk i j h (k i x)j li (k i D*)rp 123 (43) 



Performing the integrals for all the terms, with n — > —3, we 
find that there is considerable though incomplete cancellation 
between the 10 terms 



2irA 



k z 
3 re l 



- 16+ 3 



2ra+3 in 
ft 3 

3 

2 



6 + 36 + 4 + 2-18-18-12 + 18 



7,2n+3 i,n 

2wA 3 *i ^ 



(38) 



We draw attention to how the full evaluation gives a result 
for ^123 that is of opposite sign and one-sixth magnitude 
(8/3) of the value if we only consider the constant term (- 
16) for this squeezed collinear configuration. The sign of the 
mode-coupling integral is important as the bispectrum is not 
a positive-definite quantity (unlike the power spectrum) and 
the observed limits on the bispectrum may also be asymmet- 
ric about zero. Then the Case II - squeezed collinear mode- 
coupling integral is 



123 



2(4) 3 ^(n + 3) 2 V/ 



2n+3 



k G 



A- 

(39) 

As the n-dependence of 1P123 is identical for both Case I and 
Case II we can write 



123 



where 



K 



2 (4) 3 -^(n + 3) 2 



2n+3 



ka 



vl 



K = 



2 m Case I 
I Case II 



(40) 



(41) 



and we have introduced the Gaunt integral 



ym, rn n m.. 



(44) 



We substitute Eq. ( f40b into Eq. d43l for the reduced bispec- 
trum. The integrals over k 2 can be immediately done using 
the spherical Bessel function closure relation 



k 2 2 dk 2 j h {k 2 x)j h {k 2 D*) = (tt/2x 2 )5(x - D*). (45) 



and the delta function makes the x-integral trivial. We are then 
left with integrals over k\ and k% given by 



= k ( 5* 

TT 



2 (4) 3 ^ 7 (n + 3) 2 Vj. 

n+3" 



2(n+3)' 



(46) 



The fc 3 and fci integrals are similar in form to the usual Sachs- 
Wolfe term (for power-law magnetic spectra) and can be eval- 
uated analytically for any n (Eq. 6.574.2 of [66]) in terms of 
Gamma functions. For n — > — 3 we have 



K 





3 


U 2 . 


)?[ 



2 (4) 3 ^ 7 (n + 3) 2 V^. 
1 



'2Z 3 (/3 + l)2Zi(fi 



1 



(47) 



h(h + l)h(h + l)\i 2 i 3 = K TZ P 3 (|) 2 (4) 3 (n + i) 2 Vl 

(48) 
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where a factor of 1/ (D*kG) 3 ^ n+3 ^ that also appears has been 
approximated to unity for the case n —> —3 of a scale- 
invariant magnetic field index. 

For Case I - The bispectrum considering only constant an- 
gular terms 



h(h + l)h(h + l)b Wa 



1.:-sr, .in ■>■<■) 



<,^7^V(49) 



0.2 



where we have used tb ~ 10 14 GeV corresponding to a pos- 
sible inflationary epoch for magnetic field generation. This 
evaluates in different configurations to 

h(h + l)h(h + l)\l 2 l 3 « 



1 > \ 0.2 



local isosceles 

5.7 equilateral 

8.1 squeezed collinear 

9.2 midpoint collinear 

(50) 

We see that the constant-term only evaluation gives a bispec- 
trum of order 6 — 9 x 10~ 16 independent of the exact config- 
uration. The exception is the local isosceles case (where the 
constant term sum happens to vanish) but if a full evaluation 
over all angular terms is done, it yields a non-zero bispectrum. 

For Case II - The squeezed collinear bispectrum consider- 
ing all angular terms 

16 \ 



h(h + l)h(h + l)b Wa » (-1.35 x 10- 16 ) 



// +3\ 2 f ' 



0.2 



(51) 



The full evaluation of all angular terms for the squeezed 
collinear case produces two important changes (compared to 
squeezed collinear in Case I evaluation): the sign of the bis- 
pectrum has changed to negative (angular terms in mode- 
coupling contribute with different signs and change the total 
sign) and it's value had diminished by factor of 6. 

Note that the values of the reduced bispectrum given by 
Equations < f50b and ( Bit are ~ 10 6 times larger than the values 
obtained for the compensated scalar mode in SS09. This is 
essentially due to the large value of _4 3 as discussed in section 
4. 

A numerical evaluation of the mode-coupling integral in- 
volves an integrable singularity at s = and s = — k\. 
The singularity is integrable (even without a large-scale cut- 
off) provided we consider the mode-coupling integral over s 
together with the fci and k% integrals. We defer this numeri- 
cal investigation for later work and proceed with our analytic 
results. 



VI. CONSTRAINT ON PRIMORDIAL MAGNETIC FIELD 

In order to put constraints on the primordial magnetic field, 
we compare our magnetic reduced bispectrum with the re- 
duced bispectrum that arises from non-linear terms in the 



gravitational potential. These have a value given by 

h(h + lMh + l)\y 3 ~ 4 x 10- 18 /jvl (52) 

characterized by Jnl ( 15011 '). 

The reduced bispectrum for the magnetic compensated 
scalar mode lf52l 15311 was a factor of a few times 10 4 smaller 
than the standard inflationary signal (for models with Jml ~ 
1). However, we find that the magnitude of the correspond- 
ing reduced bispectrum for the magnetic passive scalar mode, 
calculated above, is ~ 30 times larger (Case II) or even up to 
~ 200 times larger (Case I) than the inflationary signal. Both 
these magnetic bispectra values are for field strength of 3 nG 
and n = —2.8, close to scale-invariant spectrum. 

Conversely, we can put upper limits on the primordial 
magnetic field by equating the magnetic bispectrum to the 
observed upper limit for the inflationary bispectrum. This is 
evaluated using the latest WMAP7 limits (95 % CL) on non- 
Gaussianity in the observed CMB -10 < < 74, 
taking the appropriate side of these limits for the different 
signs of bispectrum. Since we have considered only the 
Sachs-Wolfe contribution to magnetic field-induced temper- 
ature anisotropies, our magnetic bispectra will be accurate 
only on large scales. However, the /jvx we are using for 
comparison to data comes from WMAP data that has a 
maximum multipole of I ~ 750 and therefore the bispectra 
comparison is not exact. Yet our Bq limits are expected to be 
fairly robust as Bq depends very weakly (Bo oc /jy^J) on f^L 
as also discussed in SS09. We obtain the following upper 
limits for the magnetic field in different cases: 

Case I : For any bispectrum configuration considering only 
s-independent terms 



B < 3 nG 



(53) 



For individual configurations the magnetic field limits are: 
equilateral (2.7 nG), squeezed collinear (2.5 nG) and midpoint 
collinear (2.5 nG), all for tb ~ 10 14 GeV corresponding to 
the inflationary epoch. 

Case II : For the squeezed collinear bispectrum configura- 
tion, including all angular terms, we obtain an upper limit 



B < 2 nG 



(54) 



where this magnetic field limit value (2.4 nG) is derived using 
t b « 10 14 GeV and -10 < 

It is currently unclear what would be an appropriate epoch 
to adopt for the generation of a large-scale primordial cosmic 
magnetic field IBU 135143911 . We recall that via the logarithmic 
factor in Eq. ([H]), there is a weak dependence of Bo on tb, the 
magnetic field generation epoch. We find a 2-4 nG variation 
in the Bq upper limit when tb is varied between the inflation- 
ary and electroweak epochs tb ~ 10 14 — > 10 3 GeV, taken as 
the earliest and latest possible epochs of magnetic field gener- 
ation. However, as the electroweak transition is not expected 
to give a scale-invariant magnetic field spectrum, we caution 
that the choice of electroweak epoch may be too late to gener- 
ate a field with n — > —3. Hence, the variation in the magnetic 
field upper limit quoted above may be an overestimate. 
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VII. CONCLUSION 



We have calculated here the CMB bispectrum, on large an- 
gular scales, due to the passive scalar mode. This mode is 
sourced by the magnetic scalar anisotropic stress before neu- 
trino decoupling. The CMB bispectrum due to the passive 
scalar mode is more than two orders of magnitude larger than 
the bispectrum due to the primary inflation-induced scalar 
mode, for /jvl ~ 1 and Bq ~ 3 nG. It is also a factor of 
10 6 times greater than that previously calculated for the com- 
pensated scalar mode lf52ll53ll . This is the first calculation of 
a type of magnetic contribution to CMB bispectrum that can 
clearly dominate over the primary bispectrum at large scales. 
In the CMB power spectrum, by contrast, the passive scalar 
mode signal is highly sub-dominant to the inflationary signal 
i4(ill . Our work thus shows that the magnetically induced sig- 
nals, being intrinsically non-Gaussian, could be more easily 
distinguished from the primary CMB anisotropics, when one 
considers the bispectrum rather than the power spectrum. 

Using the WMAP7 limits for /jvx we have placed an upper 
limit on the strength of the primordial magnetic field Bo < 2 
nG which is more than an order of magnitude stronger than 



the limit from the compensated scalar mode (35 nG) obtained 
in SS09. The passive mode limit is only weakly (logarithmi- 
cally) dependent on the epoch of magnetic field generation tb ■ 
We note that for scale-invariant magnetic spectra n — > —3, the 
bispectrum and the derived magnetic field limit are indepen- 
dent of both the assumed cosmological model as well as our 
choice of smoothing scale fcg for the stochastic magnetic field 
B . 

The primordial magnetic field limit is expected to improve 
if the ISW effect and full scalar anisotropies are taken into 
account. We also expect stronger constraints to follow from 
the tensor and vector mode bispectra that dominate at low 
and high I, respectively. Better observational limits on CMB 
non-Gaussianity from PLANCK data are expected to further 
tighten primordial magnetic field limits. 
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